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Quantum entanglement of two flux qubits induced by an auxiliary SQUID 

T. Shi and Z. SorJ] 

Department of Physics, Nankai University, Tianjin 300071, China 

We revisit a theoretical scheme to create quantum entanglement of two three- levels superconduct- 
ing quantum interference devices (SQUIDs) with the help of an auxiliary SQUID. In this scenario, 
two three-levels systems are coupled to a quantized cavity field and a classical external field and 
thus form dark states. The quantum entanglement can be produced by a quantum measurement 
on the auxiliary SQUID. Our investigation emphasizes the quantum effect of the auxiliary SQUID. 
For the experimental feasibility and accessibility of the scheme, we calculate the time evolution of 
the whole system including the auxiliary SQUID. To ensure the efficiency of generating quantum 
entanglement, relations between the measurement time and dominate parameters of the system are 
analyzed according to detailed calculations. 
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Recently, superconducting quantum circuits (SQCs) based on superconducting quantum interference devices 
{Sj ' (SQUIDs) [l|, [2, [H, [H HI , have appeared to be among the most promising candidates for quantum information pro- 
cessing. Especially, schemes of quantum information transfer (QIT) and quantum entanglement by using SQC in a 
microwave cavity are carried out in many experiments E, ToL 11 . l3 |. The essential idea beyond these schemes is listed 



I. INTRODUCTION 



as follows. Firstly, one should make several SQUIDs [13l ] coupling to a single mode quantum field in a microwave cav- 
I ' ity. Secondly, QIT in the quantum computing can be realized by adjusting the external field to carry out a sequence 
j— ( , of quantum logic operations on SQUIDs. 

In this aspect, there are many experiments by using SQUIDs to realize quantum entanglement of macroscopic 

quantum devices. Here, a theoretical protocol has been proposed based on the experimentally accessible parameters 
r"| ■ [a & & Uli m which two SQUIDs are placed into a microwave cavity. Each flux qubit has a A-type energy level 
Qh configuration. An external classical field is applied to this system. Because these two three-level artificial atoms (or 
1 SQUIDs) couple to a single-mode quantum field simultaneously, an effective coupling between two three-level systems 
are induced by this single mode field. This field behaves as a data bus to link two qubits coherently. When two 
three-level systems coupling to two fields, a classical field and a quantized cavity mode, there is usually a dark state, 
q-( which contains a component with maximal entanglement of two qubits. To realize the maximal quantum entanglement 
between two qubits via the dark state, an auxiliary SQUID is used to carry out a projective quantum measurement. 
As a result of post selection measurement by wave function collapse, a quantum state with maximal entanglement 
can be obtained. 



Although this scheme is very subtle in the theoretical setup, there are still some further considerations of physical 
mechanism needed. Especially, the influence of the auxiliary SQUID should be taken into account since it should be 
treated as a quantum subsystem. According to the contemporary theory of quantum measurement, apparatus are 
usually considered as a quantum subsystem interacting with the measured system. In fact, the quantum effect of an 
auxiliary SQUID in Kansas scheme [7| is not considered in details. In this paper, we will stress the quantum effect 
of the auxiliary SQUID with respect to the whole quantum system. Our result shows that the back action of the 
auxiliary SQUID on two qubits cannot be ignored. 



II. THE GENERAL DESCRIPTION OF QUANTUM ENTANGLEMENT SCHEME OF SQUIDS 

This paper is based on the theoretical proposal of Ref . 0] ■ We first briefly review the essential idea of this paper 
with the setup being illustrated in Fig. 1. In each SQUID (I and II) the transition from the ground state |0) ; to 
the excited state \a) l resonates with the same single-mode cavity field. The transition from the first excited state |1), 
to the excited state \a} t resonates with a external classical field. The Hamiltonian H a 0, 0, EH is written in the 
interaction picture as 

H = J2 G»c|a>,<0| + Qt |o),(l| + H.c), (1) 
i=i, ii 




FIG. 1: Schematic illustration of the setup, (a) Two SQUIDs and a auxiliary SQUID are put in a microwave cavity. Here, the 
black line (dashed lines) denotes single mode cavity quantum field and grey lines denote microwave classical fields, (b) Black 
lines denote single mode cavity quantum fields. Grey lines denote microwave classical fields. Coupling constants between the 
cavity field and two SQUIDs are gi (I = 1,2). Coupling constants between the external microwave field and two SQUIDs are fi; 
(1 = 1,2). 



where gi (I — I, II) are coupling constants between the cavity field and two SQUIDs respectively, fti (I = I, II) are 
coupling constants between the external field and the two SQUIDs respectively, a and c are creation and annihilation 
operators of the cavity field mode. In Ref. Q, is time dependent to investigate the dynamics of a quantum 
information transfer between two SQUIDs. In this paper we only concern the quantum entanglement of SQUIDs and 
the effect of the auxiliary SQUID. So fi; is a constant throught this paper. 
For the Hamiltonian ([1]), the total excitation number 

Afo = (\ah(a\ + + etc, (2) 

l 

is conserved, i.e., [Hoi-A/o] = 0. The Hamiltonian can be diagonalized in each invariant subspace. In this paper, we 
focus on Afo = 1 subspace, which is spanned by the following basis 



IV>0 


= |0) 7 


|0)„|1) CJ 


IV*) 


= \a)j 


|0> w |0) o , 


IV*) 


= 10), 


l«)//|0) c , 


IV*) 


= |l)z 


|0>„|0> c , 


|V> 5 ) 


= 10), 


|l)i/|0) c , 



where |1) and |0) c denote cavity mode states with 1 and photon number respectively. Diagonalizing the matrix of 
Ho in this subspace we obtain the eigenstate 

\d) = N(n 2 g 1 \l) I \0) H \0) c + n 1 g 2 \0) I \l) H \0) c (4) 
-Q 1 Q 2 |0) J |0) Ji |l) c ) ) 




FIG. 2: A two-level system illustrating an auxiliary SQUID. The coupling constant between the cavity field and auxiliary SQUID 
is g'. 

with vanishing eigenvalue. Here, N = [(f^ffi) 2 + ) 2 + (^i^2) 2 ] -1 ^ 2 is the normalized factor. This is the so-called 
dark state because it does not couple to excited states |a),. When a photon is detected, this state is collapsed to state 
r?25i |0)jj |0) c + rii<72 |0)j 1 1) jj |0) c , which is the maximally entangled state in the case f^gi = ^1.92- This is the 
main mechanism for the generation of entanglement. In Ref. [7|, an auxiliary SQUID is employed to probe the state 
of the cavity. In this scheme the state of cavity field is measured by observing the state of the auxiliary SQUID. 

This task can be accomplished in the following process. A two-level auxiliary SQUID (see Fig. 2) can be prepared 
initially in the ground state \g) A - It is obvious that, if the cavity field is initially in a single photon state |l) c the 
auxiliary SQUID will evolve into the excited state \e) A after a half period of Rabi oscillation; if the cavity field is 
in the photon state |0) e the auxiliary SQUID always keeps in the ground state \g) A - Therefore, after a half period 
of Rabi oscillation, one can measure the auxiliary SQUID to detect the cavity indirectly. If the auxiliary SQUID 
collapses to the ground state \g) A , the cavity field collapses to the corresponding state |0) c . Then an entangled state 
of two SQUIDs 

\G) = M p 29l \l) z |0>„ + n l92 \0)j \1) H ] , (5) 

is created with M — [(SI2.91) 2 + (fii.52) 2 ] -1 ' 2 being the normalized factor. 

This the central idea of Kansas scheme, from which we can see that a post-selection measurement plays a key role 
in creating the flux qubit entanglement. The post-selection measurement is performed by an auxiliary SQUID. In the 
following, we consider the auxiliary SQUID as a part of whole quantum system. 

III. DYNAMICAL EVOLUTION 

In this section the auxiliary SQUID will be taken into account as a quantum subsystem. And the dynamical 
evolution of the measured system (two SQUIDs) will be considered. Together with the auxiliary SQUID, the whole 
system contains two parts, the auxiliary SQUID as the apparatus and the two SQUIDs as the measured system. 
Therefore, in the interaction picture, the Hamiltonian of the whole system reads as 

H = Hq + Ha, (6) 
H A = g'c\e) A (g\+H.c. 

Ha describes the interaction between the auxiliary SQUID and the cavity field with the coupling constant g' . For 
simplicity, in the following we consider the case with the parameters 

5i = .92 = .9, 

Oi = 2 = £1 (7) 

Obviously, we have 



[H,M]=0, 



(8) 
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where 



(9) 



When g' is switched off, the dark state \d) from (j4]) does not evolve since it is an eigenstate of Hq. When the auxiliary 
SQUID is employed to measure the state of the cavity field, g' is switched on, state |d) will evolve driven by the 
Hamiltonian ©. If the auxiliary SQUID is in the ground state \g) A initially, while two-SQUID system is in the state 
\d), the initial state will evolve in the invariant subspace V s with Af = 1, which is spanned by the following basis 



(10) 



Wi) = 


10), 


\0)h\1)c\9)a 


\fi) = 




|0)„|0) c | fl )^ 


\<Pz) = 


\0)j 




M = 


10), 


|0)„|0) c |e) A 


\<Ps) = 


ll>/ 


|0> w |0) c | 5 >^ 


\<Pg) = 


10), 


|i)//|o) c |3)a 



The corresponding matrix of H in the subspace V s is 



H 



( g g 

g 

9 

g' 



g' 0\ 


o o n 
ooo 



o n o o o o 

\0 !1 / 



(11) 



Considering the symmetry the whole system due to the identity of two SQUIDs (/ and II), the 6x6 matrix can be 
block-diagonalized under the new basis 



xi + 
xi + 
xi + 



\<P3», 



^7= (IV5> - l^e)), 



i 

1^4) , 
7=2 ^ 



(12) 



(|V2) + |V3», 



|V6»- 



The matrix of the Hamiltonian H can be written as 

H = H 2 ® Hi, 

where 



Ho 



1 

1 



(13) 



(14) 



and 



H d 



( V2.9 <?' \ 

V2g O 

g' 

V ft / 



(15) 



Here, g and g' are rescaled by ft: g/Cl — > 5, <?'/Sl 
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Diagonalizing two matrices, eigenvalues are obtained as 



E. 



E 



E. 



Ef 



E 



(+ 



= -1, 
= 1, 



-J=[vV-4g' 2 + ?j]i 



(16) 



where 



In the antisymmetry basis: { 
states are expressed as 



= ~E 3 , 



r 1 = l + 2g 2 + g' 2 . 



(17) 



An 



n = 1,2} and symmetry basis: { x«/ I n = 1:2,3,4}, corresponding eigen- 



N 2 ( 
N 2 ( 
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(-) 



E\- 
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9 
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9 
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1 


9 




El- 


1 


9 



x[ +) ) + V2E 3 
X ( ^)-V2E 3 



Xi +) ) + ^ 
X^)~V2E 5 



X 2 ) 

x 2 +) ) 
xi +) ) 
x 2 +) ) 



E, 


El 




9' 2 


x 3 +) 


a' 
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E 3 


El 




9' 2 


x 3 +) 


9' 




9 




E 5 


El 




9' 2 


xi +) 


9' 




9 






El 




9' 2 


xj +) 


9' 
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V2 
V2 
V2 



xi +) )) 
xi +) )) 



xi +) )) 



xi +) )) 



(18) 



where Ni and N 2 are normalized factors. 

Now, the dynamical evolution of the initial state 



|d(0)> = ^( fl |l) I |0) JJ |0) c + 9 lO^ll^lO),- 10), |0)„|l> c )| ff >. 



NiV2g 



xi +) 



(19) 



where N = 1/ \j2g 2 + 1 normalized factor, can be calculated. From (fT9|) we find that |rf(0)) is in the symmetrical 
subspace. At time t, it evolves into 



\d(t)) = F 1 (t) 



Xi +) 



F 2 (t) 



x< +) 



F 3 (t) 



xi +) 



•F 4 (t) 



xr } ) 



= F 1 (t) 10), |0)„ |1) + F 2 (t) |D) |0) o | 5 ) A 
+^3 (t) \C) |0) c | 5 ) A + F 4 (t) 10), |0)„ |0) c |e) A , 



(20) 



where 



F 1 (t) = 2iV (t) [iV 2 ^ -X cos (£ 3 *) + ' 



cos (E 5 t)], 



9 9 
F 2 (t) = -i2V2N (t) [NlE 3 X sin (Eat) + ZE 5 Y sin (E 5 t)], 
F 3 (t) = 2V2N (i) [N?X cos (E 3 t) + ZY cos (E 5 t)], 



(21) 



F A (t) = i2N(t)[-± 7 



N?E 3 El-g' 2 v ^^^ , _ 



X sin (B 8 t) + ZE 5 J * y sin (£ 5 t)], 
3 53 
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and 



X = 2g-^ . 

9 

El-1 



Y = 2g 
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Z = \ - Nl (22) 



\D) = (l«>z |0> ZJ + |0> x |a) JZ ) , 
|C) = ^(|l>zlO> JZ + |0) z |l) jrz ), 

where N (t) is the normalized factor. 

Some observations follow: (1) States \g) ^ and \e) A of the auxiliary SQUID are entangled with states of two SQUIDs 
and the cavity field; (2) When the measurement to the auxiliary SQUID is carried out, the auxiliary SQUID will 
collapse into a quantum state. If the auxiliary SQUID collapses to the ground state \g) A , one can not determine 
which state two SQUIDs will be in. This is because there exist three states (|0) 7 |0} zz , |C), and \D)) of two SQUIDs 
entangled with state \g) A - In the following section, we will discuss how to create the entangled state |C) by choosing 
optimal parameters. 



IV. CREATING MAXIMALLY ENTANGLED STATE BY OPTIMIZING PARAMETERS 

There exist two entangled states \C) and \D) in state \d(t)}. In this paper, |C) is the target state in accordance 
with Ref. The entangled state \D) contains excited states |a) ; of two SQUIDs, which couple to the environment 
with a higher probability of transitions to the first excited state or the ground state |0) z , so state \D) is unstable. 
Therefore, state \C) is a good candidate for quantum information process. 

To create entangled state |C), we first choose a special instant to to satisfy 

P (t ) = P 2 (to) = 0, (23) 

where Pi (t ) = |Pi (t )\ 2 and P 2 (to) = \F 2 (to) | 2 ■ The above equations can be regarded "entanglement condition", 
which introduce the relation between £3 and P5. The parameters can be determined when P3 (to) — I-F3 (to)| 2 , which 
is the probability of state \C) |0) c \g) A , reaches its maximum. Then the entangled state |C) can be generated with 
higher success probability. Based on the entanglement condition (|23[) , numerical simulations are employed to find the 
maximal value of P3 (P3). Numerical results are listed in Figs. 3, 4. 

As shown in Fig. 3, when Pi (to)+P 2 (to) is less than 10~ J (j = 1 ~ 6), the maximal values of P3 (to) and magnitudes 
of g and g' corresponding to these P3 (to) are given respectively. Here, the influence of the auxiliary SQUID is taken 
into account as a quantum subsystem. As a result the parameter g of the two SQUIDs and the parameter g' of the 
auxiliary SQUID must match so that the entangled state we need can be induced. The numerical results show that 
as Pi (to) + P2 (to) approaches to vanishing, the optimal area in gg' plane, within which P 3 (t n ) closes to unit. 

We also investigate the dynamical process of the entanglement creation by numerical simulation. For the case with 
fixed g and g' , under the condition Pi (to) + P 2 (to) < 1CU 6 , the time evolution of P3 (t) is calculated. 

Three groups of g and g' with (g, g') = (0.25, 1.89), (2.95, 1.10), (0.60, 1.37) are chosen to demonstrate the dynamical 



process of the creation of entanglement. Corresponding initial states are \d(0)) = N [V^g 



) ) with 



N = 1/ \/2g 2 + 1. In each case, the component of maximal entangled state in the initial state is 2g 2 /N 2 . Numerical 
results are plotted in Fig. 4. Taking Fig. 4(b) as an example, the initial state is the dark state of the two SQUIDs, 
which contains a small component (~ 0.11) of the maximal entangled state. When the coupling to the auxiliary 
SQUID is switched on, i.e., g' = — ► 1.10, to probe the state of the cavity field. After to, the quantum state which 
contains a maximal component of state \C) is generated. Then the post selection measurement to the auxiliary SQUID 
is carried out. As a result of maximizing the probability P3 (to), there exits very large probability to induce the state 
|C). This can also be found from Fig. 4(a,c). It indicates that the smaller component of state |C) is contained in 
the initial state, the longer time period to is taken. Numerical results also show that period to at which P3 (to) gets 
maximum is no longer ir/g' due to the back reaction of the auxiliary SQUID to the two SQUIDs. 



7 



p,(g.g')+p 2 (g,g')<o.i P,(g,g')+P 2 (g.g')<i°" 2 P,(g,g)+P 2 (g.g')<io 




FIG. 3: The relation between maximal values of P3 and g as well as g' when P1+P2 < W~ l (i = 1,2, ...6). The influence of the 
auxiliary SQUID is taken into account as a quantum subsystem. As a result the parameter g of two SQUIDs and the parameter 
g' of the auxiliary SQUID must match so that the target entangled state can be induced. Numerical results show that the closer 
to the entanglement condition are the results the less is the number of g and g' . 




FIG. 4: The relation between time and P3 (t). In the case of Pi (to) + P2 (to) less than 10 6 , three groups g and g' are chosen. 
Numerical results show that to which make P3 (to) maximal is not n/g' due to the back reaction of the auxiliary SQUID to two 
SQUIDs. 
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V. SUMMARY 



In this paper, the back action of the auxiliary SQUID to the two SQUIDs was taken into account in the theoretical 
scheme for creating the entangled state. The dynamical evolution of the two SQUIDs dark state ([4]) driven by the 
Hamiltonian of two SQUIDs together with the auxiliary SQUID is calculated precisely. By using the entanglement 
condition, the relation between parameters (g and g') and the evolution time to is determined. At instant to, the 
system reaches state \d(t)) — F 3 (t) \C) |0) c \g) A + F4 (t) |0) J |0) 7J |0) c \e) A . By optimizing the parameters (g and g'), 
the probability P3 (to) = I-F3 (*o) | 2 of the entangled state \C) is maximized at to. After the post selection measurement 
for the auxiliary SQUID, entangled state \C) is feasibly induced. Numerical results show that the period to given 
by the entanglement condition is not it / g' . This is because the auxiliary SQUID is taken into account as a quantum 
subsystem. 
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